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ABSTRACT 

This paper is devoted to proving that, in QCD, the lightest glueball state must be the 
scalar with J PC = ++ . The proof relies upon the positivity of the path integral measure 
in Euclidean space and the fact that interpolating fields for all spins can be bounded by 
powers of the scalar glueball operator. The problem presented by the presence of vacuum 
condensates is circumvented by considering the time evolution of the propagators and 
exploiting the positivity of the Hamiltonian. 
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In this paper I shall show that, if glueball states exist, then the lightest one must be 
the ++ scalar. There has recently been a renewed flurry of interest, both experimental and 
theoretical, in these very interesting states and the situation is beginning to clarify []I]-[|5||. 
In spite of this, the situation still remains unresolved and somewhat ambiguous so exact 
results such as that presented here are of some interest. Much detailed analysis has now 
been performed on a large amount of recent experimental data with the result that a few 
rather good candidates have emerged particularly in the region 1.5-1.7GeV 0. Potential, 
bag and instanton gas models suggest that the lowest state should be a scalar and 
that its mass should be in the above range. All of these models, in spite of having the virtue 
of incorporating the correct low energy physics of QCD, are only effective representations 
of the full theory, and so their accuracy is difficult to evaluate. However, recent lattice 
simulations of QCD based on an extensive amount of data are in general agreement with 
these model results ||. On the other hand, estimates from a field theoretic model 
indicate that the 2 ++ tensor should be the lightest state whereas a QCD sum rule analysis 
indicates that it should be the h pseudoscalar ||. This disagreement between the QCD 
sum rules and the lattice measurements is somewhat surprising since they ought to be the 
least model dependent and therefore the most reliable. However, the lattice simulations 
use a quenched, or valence, approximation, which is not generally believed to be a major 
source of error, and the QCD sum rules have difficulty satisfying a low energy theorem. In 
any case, as already stated above, the claim of this paper is that, regardless of the model 
or approximation used, the scalar must be the lightest glueball state. I shall now show 
why this must be true. 

To begin I shall first review some standard formalism as it applies to scalar and 
pseudoscalar glueballs before generalizing to arbitrary states. These spinless states can be 
described by the operators: 



G(x) = f G F^(x)Fr(x) 



and 



G(x) = f & F^(x)Fr(x) 



(1) 



where F^(x) = 



he^ 1 "*! 3 F a p(x) is the dual field tensor and fo and /<§ are constants. The 



scalar correlator 



r(x,*) = <0|T[G(x,t)G(0)]|0> 



(2) 



has a standard path integral representation: 




(3) 
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A sum over quark flavors is to be understood. By inserting a complete set of states \N) 
this can also be written 



r(x,t) = ^|(0|G(0)|JV)| 2 e^*-* PN - x ^(t) + (t -> -t) (4) 



N 



from which a corresponding Kallen-Lehmann representation can be inferred. 
A useful subsidiary quantity to consider is (for t > 0) 



Q(t) = J d 3 xT(^t) (5) 
= ^|(0|G(0)|iV)| 2 5( 3 )(p N )e lM -* (6) 



N 



where Mjv is the invariant mass of the state \N). The Euclidean version of this (effectively 
given by taking t — > ir) implies that, when r — > oo, 

Q e (t) = Q(ir) « e~ M ° T (7) 

where Mo is the mass of the lightest contributing state. An analogous result can be derived 
for r(x, t) via its Kallen-Lehmann representation where the exponential decay arises from 
the large r or \x\ behaviour of the free Feynman propagator. There are a couple of points 
worth remarking about this before proceeding. First, in pure QCD, where the ±+ glueballs 
are expected to be the lightest states in their respective channels, Mq = Mq or M G . In 
the full theory, however, the lightest states are those of 2 pions and 3 pions, respectively, 
and the glueballs become unstable resonances and mix with quark states. In that case 
Mq = Mi-k or Ms-k- On the other hand, in the limit when r becomes large, but remains 
smaller than ~ 2Mq/Fq, where is the width of the resonance, it can be shown that the 
exponential decay law, eq. ([?[) , still remains valid but with a mass Mq given by Mq rather 
than M-2k\ ( a similar result obviously also holds for the pseudoscalar case). The point is 
that, if there are well-defined narrow resonant states present in a particular channel, then 
they can be sampled by sweeping through an appropiate range of asymptotic r values 
where they dominate, since r is conjugate to Mjq @. 

The basic inequality that we shall employ is that, in Euclidean space, 

(F^±FD 2 >0 => f^G(x)>f^G(x) (8) 

Although this inequality holds for classical fields, it can be exploited in the quantized 
theory by using the path integral representation, eq. (|3|), in Euclidean space where the 
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measure is positive definite. The positivity of the measure has been skillfully used by 
Weingarten || to prove that in the quark sector the pion must be the lightest state. Here, 
when combined with the inequality (|8|), it immediately leads to the inequalities (valid for 
r > 0) 

fa 2 T E ( Xl r) > /r 2 f B (x, r) and f^Q-E^) > f~ 2 Q E (r) (9) 
By taking r large (but < 2Mq/Y'q) and using (^), the inequality 

M G < M 6 (10) 

easily follows. In pure QCD where these glueballs are isolated singularities, their widths 
vanish and the limit r — > oo can be taken without constraint. 

Although this is the result we want, its proof presumes the absence of a vacuum 
condensate E = (0|G(0)|0). It is generally believed that E ^ so the lightest state con- 
tributing to the unitarity sum in eq. ([|) is, in fact, the vacuum in which case Mq = and 
the large r behaviour of r^(x, r) is a constant, E 2 , rather than an exponential. Thus, the 
inequalities @ are trivially satisfied for asymptotic values of r since there is no condensate 
in the pseudoscalar channel. To circumvent this problem it is prudent to consider the time 
evolution of either Q(t) or T(x, t) since this removes the offending condensate contribution. 
For example, (for t > 0) 

^ = |(0|G(0)|iV)| 2 5( 3 )(p N )zM i ve lM ^ (11) 

The vacuum state clearly does not contribute to this so, in Euclidean space, the large r 
behaviour of Qe(t) is, up to a factor —Mq, just that of eq. (0). Now, (for t > 0), consider 
the following: 

^M = (0|^ G (0)|0) (12) 
= (0\i[H,G(x,t)]G(0)\Q) (13) 
= -i(0\G(x,t)HG(0)\0) (14) 

where, in the last step, the condition H\0) = has been imposed. At the classical level H is 
positive definite. We can therefore repeat our previous argument by working in Euclidean 
space and combining the inequalities (|8|) with a path integral representation for ([14]) to 
obtain (for r > 0) the inequalities 

2 dr E (x,r) 2 af £ (x,r) 2 - 2 '~ 

fa q- t > fa and f g Qe{t) > 4 Q E (V (15) 



The large r limit then leads to 



fa 2 M G e~ MGT > f^M & e~ M o- (16) 

from which (|TC)| ) follows even in the presence of condensates |J. 

Although this argument is essentially correct, it still remains incomplete in that we 
need to clarify the nature of the path integral representation for ( |i~4"D . The point is that the 
Hamiltonian, H, that generates time translations must be expressed in terms of canonical 
momentum and co-ordinate field variables. Thus, a path integral representation for (|b3| ) 
must first be written in Hamiltonian form; unfortunately, however, the resulting measure 
is not necessarily positive definite even in the Euclidean region. Thus, for our above 
argument to be valid we need to show that after integrating the Hamiltonian form over 
canonical momenta, the classical expressions for the operators can still be used and that, 
in Euclidean space, the resulting measure of the Lagrangian form remains positive. I shall 
first sketch how this comes about in a quantum mechanical context before generalizing to 
field theory. Some of the subtleties encountered are closely related to (normal) ordering 
problems that arise when dealing with operators which, like H, depend on both P, the 



canonical momentum and Q, the canonical coordinate [T(|. Consider, first, the quantum 
mechanical analog of (|j): 

r(t) = <o|r[G(t)G(o)]|o> (17) 

By analogy with eq. (|1]) the operator G is to be considered a function of Q. A path integral 
representation for this can be generated using the standard procedure of dividing up the 
infinite time interval into discrete infinitesimal sequences of size e and, at each discrete time 
t n , say, exploiting the completeness of momentum and coordinate eigenstates (labelled by 
p n and q n , respectively) fllCfl : 



r(t) = / • •/ ^/ • •/ ^e^^^-^G^ 



(18) 



Here, e = t n+ i —t n and k and I are defined such that t k =t and t\ = 0. It is understood that 
the limits N — > oo and e — > are to be taken in such a way that the total time interval 
Ne — ► zoo in order to pick out the ground state expectation value l lOfl . An analogous 
expression for dT(t)/dt can be derived in a similar fashion: 

dt J 2tt J 2tt J J 

G(q k )H(p m ,q m )G( qi ) (19) 
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Here m is defined such that t m C This freedom in the choice of t m is simply a 

reflection of the time-invariance of the operator H. If t m lies outside of this domain, then 
the integral will vanish since one can then "undo" the path integral using the completeness 
relations and move the operator H to later and later times (or earlier and earlier ones, as 
appropiate) until it annihilates on the vacuum. A similar procedure can be used to verify 
that this same expression, eq. (^9|), can be derived by directly differentiating the right- 
hand-side of eq. (p~8| ) with respect to tfc(= t). This brings down a factor i [H (pk+i, Qk+i) — 
H(pk,qk)} leading to the path integral manifestation of eq. ([Of). However, depending 
on the time ordering, one of these factors of H can be moved (again using completeness 
to "undo" the path integral) so that it eventually acts on the vacuum and vanishes. In 
this way, one can verify that the operator manipulations going from eqs. ( |T2"D - (|14| ) are 
faithfully reproduced by the path integral and that the time evolution of T(t) is given by 
eq. (0). 

Typically, and in particular in QCD, H is quadratic in P so the integrals over the p n in 
eq. ([18]) are simple gaussians which straightforwardly lead to the conventional Lagrangian 
form (up to an overall vacuum- vacuum amplitude normalization constant): 



r(t) = J d gi ...j dq N e^tl L ^), G{qk)G{qi 



(20) 



Here q' n = (q n — q n -i)/e. When elevated to field theory eq. (|2~0| ) becomes eq. (H). A 



similar set of manipulations can be carried out on eq. (|T9|). Now, however, there is an 
added complication since H occurs not only in its conventional place in the exponent of 
the measure, but also in the integrand itself. Thus, in addition to the usual gaussian 
integral, its counterpart weighted with ^p 2 n is also needed. Upon integration this leads to 
a factor (hq^ + H ) where H is a (positive) constant ~ e -2 . This background energy can 
be removed by a redefinition of H (P, Q) which only changes the energy of the vacuum, Eq, 
without affecting mass differences 0. With this definition, (|19D becomes 

dT(t) 
dt 



J dq x ...J dque'^i L ^' q - )e G(q k )H(q m ,q m )G( qi ) (21) 



As a check on the manipulations leading to this, notice that it is consistent with the 
direct differentiation of (^0|) which was normalized to the vacuum-vacuum amplitude. The 



Euclidean field theoretic version of (|2T| ) can now be used to justify the inequalities leading 
to our central result eq. (pT0|) . 



The extension of the above argument to the general case showing that the scalar 
must be lighter than all other glueball states, can now be effected. Introduce an operator, 
T^ a /3...( x )^ constructed out of a sufficiently long string of F® u (x)'s and F^ v (x)'s that it 
can, in principle, create an arbitrary physical glueball state of a given spin. Generally 
speaking a given T once constructed can, of course, create states of many different spins, 
depending on the details of exactly how it is constructed. As a simple example consider 
the fourth-rank tensor [IT 

T^ aP (x) = F^(x)F aP (x) (22) 

which creates glueball states with quantum numbers 2 ++ and ++ . Now, in Euclidean 

space, the magnitude of any component of F® u (x), or F£ v (x), is bounded by the magnitude 

i 

of [F^ v {x)F^ v {x)\ 2 . Hence, any single component of T^ uce/ 3(x) must, up to a constant, be 
bounded by G(x): 

T^ af3 (x) < fa l G(x) (23) 

This inequality is the analog of (H) and so the same line of reasoning used to exploit that 
inequality when proving ([!(]) can be used here. Following the same sequence of steps leads 
to the conclusion that Mq must be lighter than the lightest state interpolated by T^ap^x), 
from which the inequality 

M(2++) > M(0 ++ ) = Mq (24) 

follows. It is worth pointing out that the pseudoscalar analog of this operator can be 
similarly bounded thereby leading to the inequality M(2 ++ ) < M(2 '"). This argument 
can be generalized to an arbitrary T^up^^x) since, again up to some overall constants 
analogous to fa of eq. (|l|), it is bounded by some power (p) of G(x); i.e., for any of 
its components, T^q^...^) < G(x) p . Now, the operator G(x) p has the same quantum 
numbers as G{x) and so can also serve as an interpolating field for the creation of the scalar 
glueball. The same arguments used to prove that this ++ state is lighter than either the 
0" 1 or the 2 ++ can now be extended to the general case showing that it must be lighter 
than any state created by any T; in other words, the scalar glueball must indeed be the 
lightest glueball state. 

Finally, we make some brief remarks about the conditions under which the bound 
is saturated. Clearly the inequalities (|[) become equalities when F®(x) = F®(x) which 
is also the condition that minimizes the action and signals the dominance of pure non- 
perturbative instantons. In such a circumstance the ++ and 0" 1 will be degenerate. 
However, the proof of the mass inequality (|10|) only required (H) to be valid at asymptotic 
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values of \x\. Thus, the saturation of this bound actually only rests on the weaker condition 
that F be self-dual in the asymptotic region where it must vanish like a pure gauge field. 
Similarly, the saturation of the general inequality showing the scalar to be the lightest state 
occurs when all components of F® u (x) have the same functional dependence at asymptotic 
values of Although this is a stronger condition than required by the general asymptotic 
self-dual condition, it is, in fact, satisfied by the explicit single instanton solution. Thus, 
the splitting of the levels is determined by how much the asymptotic behaviour of the non- 
perturbative fields differ from those of pure instantons. This therefore suggests a picture 
in which the overall scale of glueball masses is set by non-perturbative effects driven by 
instantons (which produce the confining long-range force) but that the level splittings are 
governed by perturbative phenomena. 

This investigation was stimulated by some very enjoyable conversations in Corsica 
last Summer with Glennys Farrar, Stephan Narison and particularly Don Weingarten at 
Gluonium95. Whilst working on this problem I have further benefitted from discussions 
with Rajan Gupta and, most especially, with Tanmoy Battacharya. I would like to thank 
all of these colleagues for their helpful interactions and the DOE for its support. 
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